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We consider the indirect covering subtree problem (Kim et al., 1996). The 
input is an edge weighted tree graph along with customers located at the 
nodes. Each customer is associated with a radius and a penalty. The goal is 
to locate a tree-shaped facility such that the sum of setup and penalty cost is 
minimized. The setup cost equals the sum of edge lengths taken by the facility 
and the penalty cost is the sum of penalties of all customers whose distance to 
the facility exceeds their radius. The indirect covering subtree problem gen- 
eralizes the single maximum coverage location problem on trees where the fa- 
cility is a node rather than a subtree. Indirect covering subtree can be solved 
in 0(n log 2 n) time (Kim et al., 1996). A slightly faster algorithm for single 
maximum coverage location with a running time of 0(n log 2 nj log log n) has 
been provided (Spoerhase and Wirth, 2009). We achieve time 0(n log n) for 
indirect covering subtree thereby providing the fastest known algorithm for 
both problems. Our result implies also faster algorithms for competitive lo- 
cation problems such as (1, X)-medianoid and (l,p)-centroid on trees. We 
complement our result by a lower bound of f2(ralogn) for single maximum 
coverage location and (1, X)-medianoid on a real-number RAM model show- 
ing that our algorithm is optimal in running time. 
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1 Problem Definitions and Related Work 

We are given a tree T = (V,E) along with edge costs c: E — > M>o inducing a distance 
function d: V X V M>o. With each node u we associate a non-negative penalty ir(u). 
Let Y be a subtree of T. Then c(Y) denotes the setup cost of Y and is given by the sum 
See_E(y) c ( e ) °f ^ s e dge costs. A node is covered directly by Y if it lies in Y. If some 
node u is not covered directly it imposes the penalty ir(u) on Y. The direct covering 
subtree problem |KLTW9"6] asks for a subtree Y such that the sum of setup cost c(Y) 
and the total penalty Xm^v(Y) 7r ( n ) i s minimized. 



The indirect covering subtree problem goes one step further. A node u is said to be 
covered (indirectly) if it lies within a given distance from Y . Again, a penalty is imposed 
on Y if it does not cover u. More formally we assign to each node u some radius g(u). 
The penalty imposed onY by u is given as 



If U C V is a set of nodes then p(U,Y) := ^2 u ^u p(u,Y) is the penalty imposed on Y 
by U. The total penalty imposed on Y is given by p(Y) := p(V, Y). The indirect covering 
subtree problem [KLTW96] asks for a subtree Y of T such that the total cost c(Y)+p(Y), 
given by the sum of setup and penalty cost, is minimum among all subtrees of T. 

If we require that Y be a node rather than a subtree we obtain the single maximum 
coverage location problem [MZH83, SW09aj. It is not hard to see that single maximum 
coverage location is a special case of indirect covering subtree. (Scale all edge lengths 
and radii with a sufficiently large factor while leaving the penalties unchanged.) 

1.1 Related Work and Previous Results 

The multiple maximum coverage location problem allows the placement of a an arbitrary 
set of r nodes. On general graphs this problem is NP-hard [MZH83] while it can be 
solved in time 0(rn 2 ) on trees |Tam96| . This leads to an 0(n 2 ) algorithm for the single 
maximum coverage location problem on trees by setting r = 1. Kim et al. |KLTW96] 
provide a faster algorithm running in 0(n log 2 n). Their algorithm works even for the 
more general indirect covering subtree problem. Recently a slightly faster algorithm 
for single maximum coverage location with time 0(n log 2 nj log log n) has been reported 
|SW09a| . Finally, we remark that direct covering subtree can be solved in linear time 
1KLTW96] . 

1.2 Contribution and Outline of this Paper 

In this paper we show that indirect covering subtree can be solved in O(nlogn). This 
improves upon the previously best algorithms for this problem and single maximum cov- 
erage location on trees. Our result also implies faster algorithms for the (1, A)-medianoid 
problem and the (l,p)-centroid problem on trees. Specifically, we obtain an 0(n log??,) 
algorithm for (1, X)-medianoid and 0(n 2 log n log w(T)) and 0(n 2 log n log w(T) log D) 
algorithms for the discrete and absolute (l,p)-centroid problems on trees, respectively. 
Here, w(T) denotes the total weight of the tree and D is the maximum edge length. The 
previously best algorithms are slower by factor of 0(logn/loglogn) |SW09c| . 

Our algorithm employs the same dynamic programming framework used by Kim et al. 
|KLTW9"6] . However, we improve one of their core routines by using a more sophisticated 
technique to subdivide trees. This technique, called two-terminal subtree subdivision 
(TTST), is a simplification of the recursive coarsening strategy |SW09aj used for solving 
single maximum coverage location on a tree. The key source of our improvement is that 
we manage to avoid explicitly sorting the nodes according to their distances and radii 
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during the recursion, which has been necessary in the coarsening approach and also in 
the original algorithm of Kim et al. One further advantage of our algorithm is that it is 
a lot simpler than the recursive coarsening algorithm. 

The two-terminal subtree technique has proved successful also for other location prob- 
lems [SW10, SW09bJ. I believe that there are further problem classes where it can be 
applied. 

The paper is organized as follows. In Section [2] we briefly outline the algorithm of 
Kim et al. This is necessary, since our result relies on an improvement of a subroutine 
of that algorithm. The improved routine is then described in Section [3l In Section U] we 
provide a matching lower bound on the running time needed to solve indirect covering 
subtree. Finally, we discuss implications on related problems such as competitive location 
problems in Section [5j 

2 The Algorithm of Kim et al. 

In the sequel we will briefly describe the algorithmic approach of Kim et al. |KLTW96] 
for solving the indirect covering subtree problem. 

Let's first fix some conventions and notations. We assume that the input tree T is 
rooted at some distinguished node s. For technical reasons we shall adopt the convention 
that s is the father of itself. Let v be an arbitrary node. Then f(v) denotes the father 
of v. We write T v for the subtree of T hanging from v and T+ for the union of T v with 
the edge (v, f(v)). 

Kim et al. reduce the solution of the problem to the computation of the values 
p(v),p(T v ,v) and p(T v , f(v)) for all nodes v. They show that one can determine an 
optimum to the subtree location problem in linear time once these values have been 
precomputed for all nodes v. 

To convince ourselves, assume that we have computed the values p(v), p(T v ,v) and 
p(T v , f(y)) for all v £V. Then define 

C{v) := min{ c(Y) + p(T v , Y) \ Y is subtree of T v containing v } , 

and 

C + (v) := min{ c(Y) + p(T v ,Y) \ Y is subtree of containing f(v) } . 

It is not hard to see that the optimum cost can now be expressed by 

min(C(i>) +p(v) —p(T v ,v) ). 

Moreover, the C(-)- and C + (-)-values can be computed in linear time by means of a 
simple bottom-up dynamic programming approach. To this end assume that v is a leaf 
of T then 

C(v)=0 and C + {v ) = mm{p(v, f(v )), c(v, f(v))} . 
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Otherwise, we have 

c( V )= c+ ( u )< 

u is son of v 

and 

C + {v) = mm{C(v) + c(v, f(v)),p(T v ,f(v))} . 

From this equations it follows that an optimal solution can be determined in linear 
time in a bottom-up fashion once the values p(v),p(T v ,v) and p(T v ,f(v)) have been 
computed for all v € V. Kim et al. suggest an algorithm with running time 0(n log 2 n) 
to compute these values. 

This algorithm is based on the so-called bitree model. In this model, each (undirected) 
edge (u,v) of the input tree is replaced with two anti-parallel, directed arcs (u, v), (v,u). 
We call the resulting tree T' bitree of T. With each arc (u, v) of the bitree we associate 
a cost Ct'(u,v) representing the length of this arc. But in contrast to the edges of the 
input tree T we allow these costs to be negative and asymmetric. This induces a distance 
function dx> ■ V x V — > Q where dx'(u, v) is the length of the unique u-v-path in T'. Now 
we define the penalty cost p'(u,v) imposed on v by u to be zero if dT'(u,v) < g(u) and 
ir(u) otherwise. We set p'(v) = Ylu&v p' i u -> v )- 

The algorithm of Kim et al. is based on a subroutine for efficiently computing p'(v) 
for all nodes v on a given bitree T' . By means of such a subroutine it is then possible to 
calculate the values p(v), p(T v ,v) and p(T v , f(v)) for all v in the input tree T. It follows 
from the above discussion that the knowledge of these values enables us to identify an 
optimal tree-shaped facility. 

It remains to explain how we can employ such a subroutine to determine p(v), p(T v ,v) 
and p(T v , f(v)) for all nodes v of the input tree T which, in turn, is sufficient to build an 
optimal tree-shaped facility. 

First we describe how we can determine p(-). For this purpose we simply set ct> (u, v) := 
Ct'(v,u) := ct(u,v) for all edges (n, v) of the input tree T. It is then immediately clear 
that p[v) = p'(v) for all v £V. 

In order to compute p(T v ,v) for all v E V we set Ct'(u, f(u)) := ct(u, f(u)) and 
CT'(f(u),u) := — oo for all u ^ s. This construction ensures that the penalty cost 
p'(u,v) is always zero if u is not a descendant of v. Thus p(T v ,v) = p'(v) holds for this 
construction. 

Finally, we wish to determine p(T v ,f(v)). To this end we introduce on each edge 
(v , f(v)) of T a new node f'(v) such that edge (v , f'(v)) has length ct(v, f(v)) and edge 
(f'(v),f(v)) has length zero. This increases the number of nodes to 2n — 1. We set 
7r (/'(?;)) and g(f'(v)) to zero. It is easy to see that p(T v , f(v)) in the original tree equals 
p(T u , u) in the newly constructed tree where u := f'{v). Hence the problem of computing 
p(T v ,f(v)) for all nodes v can be reduced to the problem of computing p(T v ,v), which 
has been described before. 

Kim et al. provide a subroutine to compute the p'(-)-values on a bitree with n' nodes 
in 0(n' dog 2 n') time which yields immediately. 
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Theorem 1 QKLTW96J) The indirect covering subtree problem on a tree can be 
solved in 0(n log 2 n). □ 



3 An 0(n log n) Algorithm 

In this section we describe an algorithm for the indirect covering subtree problem with 
running time 0(n log n). 

Our algorithm uses the algorithmic framework of Kim et al. described in Section [2j 
Specifically, we will provide an improved routine for computing the values p'(-) on a 
given bitree in 0{n log n) which can then be extended to an algorithm with the same 
asymptotic running time for solving indirect covering subtree. 

The basic approach of the routine of Kim et al. for computing p'(-) is divide-and- 
conquer. It partitions the node set V into two sets Vi, V2 of bounded size such that both 
induce subtrees and have exactly one node (called centroid |KA75| ) in common. Then 
it sorts the sets Vi and computes, by means of a clever merge-and-scan procedure, for 
all v G Vi the penalties p'(v,Vj) of the users in Vj where j ^ i. Applying the routine 
recursively to the sub-bitrees induced by V\,V2 one can determine the p'(v,Vi) -values 
also for each v G V{. Finally, one obtains the total penalty p'(v, V) of any node v (zV hy 
adding p'(v, V\) and p'(v, V2). 

Our routine proceeds in a similar way but uses a more sophisticated subdivision, which 
allows us to avoid the explicit sorting thereby supressing the additional log-factor. Spo- 
erhase and Wirth |SW10| ISW09b| used an analogous subdivision technique for solving 
competitive location problems on undirected trees. 

Consider the (undirected) input tree T = (V, E). We may assume that T has maximum 
degree three. Otherwise, we can split nodes of larger degree by introducing suitable zero- 
length edges and zero-weighted nodes. Let T' be the bitree corresponding to T. 

If s and t are distinct nodes then T' st denotes the maximal sub-bitree of T' having s and 
t as leaves. Let V s t be the node set of T' st . We call s and t terminals and T' st two-terminal 
sub-bitree (TTSB). 

Our algorithm divides the input bitree recursively into TTSBs. Since we are dealing 
with a degree-bounded bitree we can subdivide any TTSB S into at most five TTSBs, 
called child TTSBs. Each of these child TTSBs has at most t^S*! + 1 nodes. 

Lemma 2 Let S be a TTSB with maximum degree three. Then S can be partitioned 
into at most five edge-disjoint TTSBs each of which having at most \\S\-\-l nodes. This 
subdivision can be computed in 0{\S\) time. 

Proof. Let S be a TTSB with maximum degree three and terminals u and v. Let m be 
the unweighted median of S, which can be computed in OdS"!) by means of Goldman's 
algorithm |Gol71| . (All node and arc weights are temporarily set to one throughout this 
proof.) It is a well-known fact that m has the following property: Each of the connected 
components of S — m has at most \\S\ nodes. Hence, if m lies on path P(u,v) then 
S — m contains at most three components that form the desired subdivision (confer left 
part of Figure [1]). If m does not lie on P(u,v) then consider the node m! on P(u,v) 
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that is closest to m (confer right part of Figured]). Then S — {m,m'} has at most five 
connected components. All of the child TTSBs obtained this way have clearly at most 
nodes. □ 




Figure 1: The two cases in the subdivision of a TTSB. 

Consider a TTSB T' st . We introduce the lists L^ s {T' st ) and L gs {T' s ^). Both lists contain 
all nodes v of T' st sorted in increasing order with respect to the values d,T>{s,v) and 
q(v) — dj"{v, s), respectively. The lists Ldt(Tg t ) and L e j{T' st ) are defined symmetrically. 

The algorithm computes p'{v,T' st ) for all v 6 T' st as well as the four lists L^ s (T!, t ), 
L(£t(Tg t ), L ss (Tg t ) and L Sy t{T' st ) for any TTSB T' st occurring during the recursion. We 
shall see that these information can be propagated inductively from child towards parent 
TTSBs such that we will have computed p'(-,T') = p'(-) at the top of the recursion. 

To this end consider an arbitrary TTSB S = T' st being subdivided into at most five 
child TTSBs Si with terminals Si, t{. Moreover assume that we have already computed 
p'(-) and the four lists corresponding to Si for all Si. 

We start with computing L^ S (S). To this end we maintain a list L which is initialized 
with an empty list. Now we perform the following operations for all child TTSBs Si'. 
Assume that Si is the terminal of Si closest to s. Then the list Ld tSi (Si) contains all nodes 
v £ Si with associated sorting keys dx'(si,v). Now we create a copy V of this list and 
add the value d,T'(s, Si) to all sorting keys which does not affect its order. As a result L' 
contains all nodes v of Si sorted with respect to their distance d,T'(s,v) from terminal s. 
Finally we merge L with L' . After having carried this out for all child TTSBs Si the list L 
equals the list Ld s (S) we are looking for. The list L g ^ s (S) is computed very similarly with 
the difference that we subtract the value dT'(si,s) from the sorting keys when handling 
the list Lg jS (Si). The respective lists for terminal t are computed symmetrically. The 
total running time for computing the four lists associated with S is 0(|iS|) since we handle 
a constant number of child TTSBs. 

We are now going to explain how p'(v, S) can be determined for all v 6 S. To this end 
assume that v is contained in some Si. Since we already know p'(v, Si) by the inductive 
hypothesis it suffices to determine p'(v,Sj) for all Sj ^ Si and to add these values to 
p'(v,Si). Consider an arbitrary Sj ^ Si and assume that Sj, Sj are the terminals of these 
TTSBs closest to each other. We create a copy L' of list L Q ^ Sj (Sj) and substract the 
distance dT'(sj, Si) from all sorting keys in this list. As a result L' contains all nodes u of 
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Sj sorted with respect to the key g(u) — dx'(u, Si). At this point we can compute p'(v, Sj) 
for all v £ Si by using the merge-and-scan procedure of Kim et al. To this end we merge 
the sorted list L' with the sorted list Si {Si) an d store the result in V . We assume that 
the nodes in L' are sorted in increasing order with respect to their numerical sorting keys. 
Ties are broken in favor of nodes in Si. This can be achieved in linear time 0(|<Sj| + \ Sj\). 
Now recall that a node u £ Sj imposes a penalty ir(u) on v £ Si if dx'(u,v) > g(u) or 
equivalently dj<'(si,v) > g(u) — dnpi{u,Si). This is tantamount to that u precedes v in 
V . Hence, in order to compute p'(v,Sj) for all v £ Si it suffices to traverse V once. In 
doing so, one can maintain the sum of penalties of all nodes u £ Sj encountered so far, 
which equals the penalty p'(v, Sj) whenever a node v £ Si is reached. 

The running time of this merge-and-scan operation is 0(|5i| + \Sj\) since the necessary 
sorted lists have already been computed. Thus we can compute p'(v,S) for all v £ 5 
in total time Od^l) once we know the p'-values and respective lists for all child TTSBs 
of S. 

Note that the bottom of the recursion, that is, when T' st consists merely of the pair 
(s, t) and (i, s) of anti-parallel arcs can trivially be handled constant time. 

To sum up, this leads us to an algorithm whose running time /i( | S' | ) can be described 
by the following recurrence 

k 

h(\S\) = 0(\S\) + Y,H\Si\), 
i=i 

where k < 5, Xa=i 1^1 = 1*51+4 and \Si\ < ^ | ^S* | + 1 . This implies that h(n) is O(nlogn). 

Theorem 3 The indirect covering subtree problem and hence also the single maximum 
coverage location can be solved in time O(nlogn). □ 

4 A Matching Lower Bound 

In this section we complement our algorithm with a lower bound f2(n log n) on the running 
time for solving single maximum coverage location on a tree. This shows that (for certain 
computational models) our algorithm is optimal. 

We make use of a recent result which is summarized in the following theorem. 

Theorem 4 ([B AGOl] ) Let W C W 1 . If W is recognized in time t(n) on a real- 
number RAM that supports direct assignments, memory access, flow control, and arith- 
metic instructions {+,— , x,/} then t(n) = £l(\og(3(W°)). □ 

Here, W° denotes the interior of W and (3(W) denotes the number of connected com- 
ponents of some set W' C ]R n . 

To prove our lower bound we introduce a variant of the set disjointness problem. To 
this end let n £ N. The set W n C R^ 1 contains all tuples (xi, . . . , x n , yi, . . . , y n ) such 
that x\ < . . . < x n and Xi ^ yj for all pairs i,j. Consider a permutation ir on the set 
{1, . . . , n} and some tuple x\ < y^m < %2 < V-k(2) < ■ ■ ■ < x n < y^i n ) in W n . It is easy 
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to see that for different permutations such tuples lie in different connected components 
of W° so W° contains at least n\ connected components. Hence any RAM of the above 
described type takes time fi(nlogn) to recognize W n . 

We establish a linear time reduction from the problem to recognize W n to the single 
maximum coverage location problem on a tree with 0(n) nodes. To this end consider a 
tuple (xi, . . . , x n , yi, ... , y n ) for which we want to decide whether or not it is contained 
in W n . 

First we check if x\ < ... < x n . Then we create an edge (u, v) of some length c(u, v) > 
max{xj,yj | i = 1, ... ,n} and choose some radius q such that q > c(u,v). For any yi 
we create two edges (u, Ui) and (v, Vi) of lengths Q — yi and yi + Q — c(u, v), respectively. 
Finally, we create for each X{ a node X{ on edge (u,v) with distance d(u,Xi) := x^. For 
each node z in the node set V := { Ui, Vi, Xi \ i = 1, . . . , n } U {u, v } we set ir(z) := 1 and 
q(z) := q, which completes the reduction. 

First suppose that we locate a facility outside the path P(u, v). Assume that the facility 
is located at some node m. Then the distance of U{ to u is positive and d(u,Vj) > q for 
any j. Hence, none of the nodes Vj is covered by U{ and the penalty cost imposed on Ui 
must be at least n. The case where is the facility is placed at some node v i is treated 
analogously. 

Now suppose for a moment that we can locate a facility anywhere at the path P(u, v), 
that is, also at interior points of edges on P(u,v). The point x where the facility is 
located can then be identified with the distance d(u,x). First, all nodes on P(u,v) are 
covered by x since q > d(u,v). Due to our construction x covers all nodes Ui where 
x < yi and all nodes v j where x > yj. Thus, the penalty imposed on x is exactly n if x is 
not contained in the set {y\, . . . , y n }. If x = yj then x covers both Uj and Vj and hence 
the penalty is bounded by n — 1. Since the facility can only be placed at nodes Xj, that 
is, at distances xi from u we conclude that the minimum penalty cost is n if the input 
tuple (xi, . . . , x n , j/i, ... , y n ) lies in W n and n — 1 otherwise. 

Theorem 5 Any real-number RAM that complies with Theorem^takes at least fi(n log n) 
time to solve single maximum coverage location on a tree even for unit penalties and uni- 
form radii. □ 

5 Implications for Related Problems 

The variant of single maximum coverage location where the facility can be placed not 
only at the nodes but also at interior points of edges is called the absolute maximum 
coverage location problem. Kim et al. show [KLTW96] that a set of 0(n) critical points 
(that is a set of point which is guaranteed to to contain an optimal point) for absolute 
single maximum coverage location can be found in time O(nlogn). We infer that also 
the absolute variant can be solved in 0(n log n) on a tree. 

Another implication of our result leads us to the realm of competitive location. Let a 
graph G = (V, E) and r,p < n be given. We assume that the graph is edge and node 
weighted. Let X, Y C G be sets of nodes or interior points of edges. Then w(Y -< X) 
denotes the total weight ^2{w(u) \ u € V and d(u,Y) < d(u,X)} of nodes that are 
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closer to Y than to X. Given some point set X the goal of the (r, X)-medianoid problem 
|Hak83| is to identify a set Y of r points such that w(Y -< X) is maximized. This 
maximum weight is denoted by w r (X). The goal of the (r,p)-centroid problem is to find 
a p-element point set X such that w r (X) is minimized. 

By setting g(u) := d(v, X) — e (where e is a suitably small constant) and tt(u) := w(u) 
one can easily verify that (r, X)-medianoid is a special case of the multiple maximum 
coverage location problem with r servers. On general graphs the problem is NP-hard 
|Hak83| . It can be solved efficiently in 0(rn 2 ) on trees |Tam96 ] . Our result leads to 
an 0{n log n) algorithm for the absolute and the discrete version of (1, X)-medianoid on 
trees. 

Corollary 6 The discrete and the absolute (l,X)-medianoid problem can be solved in 
O(nlogn) on trees. □ 

It is not hard to extend the lower bound provided by Theorem [5] to (1, X)-medianoid. 
Since the radii of the tree constructed in the reduction are uniform and hence all equal 
some number g, we can furnish each node z on that tree with a pendant leaf z' at a 
distance d(z, z') = g + £. The set X contains exactly those pendant leaves. It is clear 
that for any node y in this enhanced tree T' the gain w(y -< X) equals exactly w(T)—p(y) 
in the original tree T. This implies that both instances lead to the same optimum. Thus 
also the algorithm for (1, X)-medianoid is optimal in terms of the running time. 

Now let's turn our view to the (r,p)-centroid problem. The problem is known to be Un- 
complete on general graphs [NSW07] and NP-hard even on path graphs |SW09c| . H ow- 
ever, both the absolute and the discrete variant of (l,p)-centroid on trees can be solved in 
polynomial time O (n 2 log 2 n log w (T) / log log n) and O (n 2 log 2 n log w (T) log D / log log n) , 
respectively |SW09c| . Here D is the maximum edge length of the input tree T. Those 
algorithms rely on 0(n log w(T)) (resp. 0(n log w(T) log D)) calls to a subroutine solving 
(1, X)-medianoid on a tree. 

The algorithm provided here allows us to solve (1, X)-medianoid in 0(nlogn) which 
yields. 

Corollary 7 The discrete and the absolute (l,p)-centroid problem for trees can be 
solved in 0(n 2 log n log w(T)) and 0(n 2 log n log w(T) log D), respectively. □ 

6 Concluding Remarks 

We have provided an 0(n log n) algorithm for solving the indirect covering subtree prob- 
lem which improves upon the previously best algorithms for this problem and single 
maximum coverage location on trees. We have also shown that our algorithm is optimal 
for certain unit-cost RAM models. Our result leads also to an optimal algorithm for 
(1, X)-medianoid and faster algorithms for (l,p)-centroid on trees. 

It would be interesting to identify larger problem classes of location problems on trees 
where the the two-terminal subtree technique can be applied. It would also be worth 
investigating the existence of faster algorithms on path graphs. 
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